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1. Introduction 

To warrant high reliability of key components of systems, stand-by units are used. 
If such a component fails and at least one not failed stand-by unit remains, it 
operates instead of the key component. If this one fails, the second stand-by unit 
is used and so on until the last stand-by unit fails. We suppose that switching is 
instantaneous and there are no repairs. 

If the stand-by units are functioning in the same "hot" conditions as the main 
unit then usually after switching the reliability of the stand-by units does not 
change. But "hot" redundancy has disadvantages because any one of the stand-by 
units fails earlier than the main one with the probability 0.5. 

If the stand-by units are not operating until the failure of the main unit ("cold" 
reserving), it is possible that during and after switching the failure rate increases 
because the stand-by unit is not "warmed" enough [1]. So "warm" reserving is 
sometimes used [2]: stand by units function under lower stress than the main one. 
In such a case the probability of the failure of the stand-by unit is smaller than 
that of the main unit and it is also possible that switching is fluent, i.e. switching 
from "warm" to "hot" conditions does not do any damage to units. What does it 
mean mathematically? 

Let us consider a system of m units: one main unit and m — 1 stand-by units. 
We shall use notation 5(1, m — 1) for such systems. 

Denote by T\, F\ and f\ the failure time, the c.d.f. and the probability density 
function of the main unit. The failure times of the stand-by units are denoted by 
T2 , . . . , T m . In " hot" conditions their distribution functions are also F\ . In " warm" 
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conditions the c.d.f. of Tj is F 2 and the p.d.f is f 2 , i = 2, . . . , m. If a stand-by unit 
is switched from "warm" to "hot" conditions, its c.d.f. is different from F\ and F 2 . 

In 1966 N.M.Sedyakin [3] formulated the "physical principle in reliability". The 
idea is the following. For two groups of units from identical populations functioning 
under different constant stresses x\ and x 2 , two times t\ and t 2 are equivalent in 
the sense that the probabilities of survival until these moments are equal: 

1 - Fi(ti) = P(Ti > *i|xi) = P(T 2 > t 2 \x 2 ) = 1 - F 2 (t 2 ). (1) 

If the first group of units is tested under the constant stress x\, the second group 
is tested under the constant stress x 2 up to time t 2 and under the stress x\ for 
t > t 2 , i.e. under the step-stress 



x\u) 



x 2 , < u < t 2 , 
Xi, u>t 2 , 



then Sedyakin's hypothesis is the following: for all s > 

\(t 2 + s\x(-)) = A(ti + s|xi); (2) 

here \(t\x) denotes the hazard function under the stress x and t\ satisfies (1). The 
equality implies 

P(T 2 >t 2 + s\T 2 >t 2 ,x(-)) =P(Ti >t 1 + s\T 1 >h, Xl ), 

which means that the probability not to fail time s under the same stress after two 
equivalent times t\ and t 2 are the same for two groups of units. So switching from 
the stress x\ to the stress x 2 does not do any damage to units. If the stresses x 2 
and x\ mean "warm" and "hot" conditions, respectively, and the switch on time t 2 
would be non-random then we could use formula (2) for computing the conditional 
distribution of the switched on stand-by unit (see [6] , [7] ) . In the case of redundant 
systems the formula (2) should be modified because the switch on times are random 
- stand-by units are switched on after the failure of the unit operating in "hot" 
conditions. 

The failure time of the system S(l,m- 1) is = T x V T 2 V . . . V T m . Using the 
principle of Sedyakin, we consider the following model for the c.d.f. of the random 
variable T. As = (Ti V T 2 V . . . V T m _i) V T m , we can consider this system as 
a system 5(1, 1) with one main element (which itself is a system 5(1, m — 2)) and 
one stand-by element. 

Denote by Kj and kj the c.d.f. and the p.d.f. of respectively, (j = 2, . . . , m), 
K\ = Fi, k\ = j\. The c.d.f Kj can be written in terms of the c.d.f Kj-i and F±: 

Kj{t) = P(T^) <t) = P^'- 1 ) < t,Tj < t) = f P(Tj < t\T^-V = y )dK j ^ l {y). 

Jo 

(3) 

We use the Sedyakin's principle modelling the conditional distribution P(7j < 
t\T^~^ = y). If T^ -1 ) = y and the stand-by unit has not failed until this moment 
then it is switched on. 
Hypothesis Hq: 



(4) 
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"equivalent to y moment" g(y) is found from the equation F\(g(y)) = F 2 (y), so 

g(y) = F 1 - 1 (F 2 (y)). 

The formula (4) implies 



Kj(t) = f F l (t + g(y) - y)dK j _ 1 (y). 
Jo 



(5) 



So the distribution function of the system with m — 1 stand-by units is defined 
recurrently using formula (5) (j = 2, . . . , m). 

In particular, if we suppose that the distribution of units functioning in "warm" 
and "hot" conditions differ only in scale, i.e. 

F 2 (t) = Fi(rt), (6) 

for all £ > and some r > 0, then g(y) = ry. 
Hypothesis Hq\ 

^-^-t/i^ + ry-yjif t>y., ^ 

Conditionally (given T^ -1 ^ = y) this hypothesis corresponds to the accelerated 
failure time (AFT) model [4], [5], [8]. Under model (7) the distribution function of 
the system is obtained using recurrent formulas 

Kj{t) = f F 1 (t + ry-y)dK j _ 1 (y). (8) 
J o 



2. Nonparametric estimation of the reliability of the redundant system 

Suppose that the hypothesis Hq is true and the following data are available : 

a) complete ordered sample Tn, . . . , T± ni of the failure times of units tested in 
"hot" conditions; 

b) the time to obtain complete data in "warm" conditions may be long, so we 
suppose that n 2 units are tested up to time t\ in "warm" conditions and the ordered 
first failure times T 2 i, . . . , T 2m2 are obtained. 

Set 

Hi n 2 
i=l i=l 



ni ?i2 

Y 1 (t) = J2 l{r 14 >t}» Y 2 (t) = l{r M > t) t< tl }- 

i=l i=l 



Note that the random variables Tu / r and T 2 i can be interpreted as order statistics 
from samples of size n\ and n 2 , respectively, from the population having the c.d.f 
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F 2 . So if we denote 



TLX 



N x (t) = l{T 1( /r<t} = Ni(rt), N 2 (t) = N 2 (t), 



i=i 



= J] l {Tl( /r>t} = Yi(rt), Y 2 (t) = Y 2 (t), 

i=l 

then Nelson- Aaalen type estimator (still depending on r) of the cumulative hazard 
function A2 = — In S 2 can be considered: 

~ [* dNiju) + dN 2 {u) [ f dNijru) + dN 2 (u) 

2[,r) Jo Yi(u)+Y 2 (u) Jo yi(ru)+y 2 (u) ' [) 

Taking into consideration that the difference 



M 2 (t) = N 2 (t) - [ Y 2 (u)dA 2 (u) 
Jo 



is a martingale on [0, t±] with respect to the filtration generated by the data, and 
~EM 2 {ti) = 0, the parameter r can be estimated using the estimating function 

U(r) = N 2 {h) - t Y 2 (u)dA 2 (u,r) = 







Y^vMdN^v) Y 2 (u)dN 2 (u) 
2(1) J Y^+Y^v/r) J Y 1 (ru)+Y 2 (u)- 

U(r) is a non-increasing step function, 

so the parameter r is estimated by the statistic 

f = U^iO) = sup{r : U(r) > 0}. 
The estimator of the cumulative hazards Ai and A 2 are 

diVi(n) , dN 2 (u) 



A l(t ) _ A^/r.r) _ J y^r^ + ) + 



, , t Y 1 (T lt ) + y 2 (T lt /f) T ^ /f Y 1 (fT 2t ) + ^(T 2l ) ' 



A 2 (t) = Ai(ft). 
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The estimators of the c.d.f. Si = 1 — Fi are the product integrals of the estimators 
Ai(i) and A 2 (i), so 



#i(t) = l-7T < a <t(l-cLii(a)) 



n i-^r-w^d n 1 



r , v Y 1 (T li ) + Y 2 (T li /r)J T ^ /f V n(rr 2i )+y 2 (T 2i )y ' 



F 2 (t) = A(rt). 

Mixing all moments Tij and fT 2 j and ordering them, we obtain the sequence of 
random variables T\ < ... < T ni+iri2 . The estimators Fi(t) and F 2 (£) can be 
written: 

The c.d.f. -ftT m of the redundant system is estimated using the following recurrent 
equations (j = 2, . . . , m): 



A(* A _ l( t, + E (10) 



rt<T,<< 

The estimator of the mean failure time of the system is 

poo nt+m.2 

ft = / tdK m (t) = Vt ; [£ m (T0 - K m (Tj_i)]. 
J ° i= i 

Failure times from exponential distribution were simulated: 

Tij ~ £/(Ai), T 2 j - £(A 2 ), A 2 = rX 1 . 

The graphs of the estimators of the c.d.f. F\ and K rn (m = 2,3,4), in the case 
of complete samples and different sample sizes are presented in Figures 1 and 2. 
Increasing the number of stand-by units increases the reliability of the redundant 
system. 



3. Goodness-of-fit 

The given estimators K m of the c.d.f. of the redundant system S(l, m — 1) can be 
used if the hypothesis Hq is true. If switching from "warm" to "hot" conditions 
does not damage units in the system £(1,1) then it is natural that this is true 
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Sample size is 50 




Figure 2. Estimators of the c.d.f. Fi and K m (m = 2,3,4), n± = n2 = 100, ti = oo 



for the system S(l,m — 1), m > 2. So we do not need a test with 5(1, m — 1) 
systems and it is sufficient to give tests for the hypothesis Hq and Hq when only 
one stand-by unit is used. 

Suppose that the following data are available : 

a) the failure times Tn, . . . , T\ ni of n\ units tested in "hot" conditions; 

b) the failure times T21, . . . ,T2 n2 of 712 units tested in "warm" conditions; 

c) the failure times T\ , . . . , T n of n redundant systems (with " warm" stand-by 
units). 

The tests are based on the difference of two estimators of the c.d.f. F. The first 
estimator is the empirical distribution function 



i=l 



(11) 
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The second is based on formula (5), i.e. 

F&(t)= f F 1 (t + g(y)-y)dF 1 (y), 
Jo 

where (hypothesis Hq) 

g(y)=F 1 -\F 2 (y)), Fj (t) = -Y^MT^t}, F 1 -\y) = mf{s:F 1 (s)>y}, 

(12) 



3 i=i 



or (hypothesis Hq) 



2 3 

g(y) = fy, r = -±, fa = — T^. (13) 

rij , =1 



The test is based on the statistic 



/•oo 

X = ^/n~ [F^\t) - P( 2 \t)]dt. (14) 
Jo 

It is a natural generalization of Student's t-test for comparing the means of two 
populations. Indeed, the mean failure time of the system with c.d.f. F is 

poo 

ix= \ [1-F(s)]ds, 
Jo 

so the statistic (12) is the normed difference of two estimators (the second is not 
the empirical mean) of the mean ji. Student's t-test is based on the difference of 
empirical means of two populations. 

It will be shown that in the case of both hypothesis Hq and Hq the limit distri- 
bution (as rii/n — > l; L £ (0, 1), n — ► oo) of the statistic X is normal with zero mean 
and finite variance cr 2 , (see Theorems 1 and 2). 

The test statistic is 



± n 



x ,2 



a 



where a is a consistent estimator of a. The distribution of the statistic Y n , is 
approximated by the standard normal distribution and the hypothesis Hq (or Hq) 
is rejected with approximative significance value a,0 < a < 0.5, if Y 2 > Xi- Q (l)> 
where Xi-a(l) i s the (1— a)-quantile of the chi-squared distribution with one degree 
of freedom, see, for example, Greenwood and Nikulin (1996). 
Let us find the asymptotic distribution of the statistic (14). 

Theorem 3.1 : Suppose that rii/n — ► Zj G (0, 1) ; n — > oo and the densities fi(x), 
i = 1,2 are continuous and positive on (0,oo). Then under Hq the statistic (14) 
converges in distribution to the normal law N(0,a 2 ), where 

i 2 2 

a 2 = VarCTO + -Var(i7(T ll )) + ^Var(T 2l ), (15) 
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where 

H(x) =x[c + r-l- F 1 (x/r) - rF 2 {x)\ + rE(l {Ti ^ x/r} T u ) + rE(l {T2i < x} T 2i ), 



1 f°° 

c=~ y[l-F 2 (y)]dF 1 (y). 

^2 Jo 

Proof : The limit distribution of the empirical distribution functions is well known: 

y/H(Fi - Ft) Z Ui, V^{F {1) -F)%U (16) 

on D[0, oo), where — > means weak convergence, U±,U 2 and U are independent 
Gaussian martingales with Ui(0) = U(0) = and the covariances 

cov(C/i(si), I7 i (a 2 )) = r Fi( Sl A a 2 )[l - ^(si V a 2 )], 



cov([/( Sl ), I7(s 2 )) = F(ai A s 2 )[l - F( Sl V s 2 )]. 

Under hypothesis Hq the difference of the two estimators of the distribution func- 
tion F can be written as follows: 

F^(t)-F^(t) = F^(t)-F(t)- f F 1 (t+g(y)-y)dF 1 (y)+ f F 1 (t+g(y)-y)x 

Jo Jo 



dF^y) = F^(t) - F(t) - /Vi(* + 9(V) -y)-F 1 (t + g(y) - y)]dF 1 (y)- 

Jo 



I [(F 1 (t + g(y)-y)-F 1 (t + g(y)-y))-(F 1 (t + g(y)-y)-F 1 (t + g(y)-y))]dF 1 (y)- 
Jo 



[ [A (t + g {y)-y)- F 1 (t + g(y) - y)\ [dF-y (y) - dF 1 (y)} - 
Jo 



/W + g(y) -y)- F 1 (t + g(y) - y)]dF l {y)- 
Jo 



f [A (t + g(y)-y)- F x (t + g(y) - y)\ [dF-y (y) - dFy (y)} - 
Jo 



1 Fy(t + g(y) - y)[dFi(y) - 



1 
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The statistic (14) can be written 



pOO 

X= / ^[F {1 \t) - F(t)]dt- 
Jo 



poo pt 

/ dt y/E[F 1 (t + g(y)-y)-F 1 (t + g{y)-y)]dF 1 (y)- 
Jo Jo 



poo pt 

/ dt y /fi[F 1 (t + g{y)-y)-F 1 (t + g(y)-y)]dF 1 (y)- 
Jo Jo 



poo pt 

/ dt F 1 (t + g(y)-y)d{MF 1 (y)-F 1 (y)}} + op(l) 
Jo Jo 



Set a 2 = Var(Tjj), j = 1,2. The convergence 



/-oo 

V^G&i - Mi) ■ E >Yj = - U 3 {y)dy ~ iV(0, of 

■/ o 

implies 



^(r - r) £ y = -(11 - rY 2 ) ~ AT(0, %f + f )). 

M2 M 2 <l £2 



Formulas (16)- (18) imply 

Jo Jo 



poo pt 

/ dt / ^[Fi (t + g(y) -y)-F 1 (t + g(y) - y)\ dF l (y) Z 
Jo Jo 

poo poo 

cYi - rcY 2 = -c / Ui(y)dy + rc U 2 (y)dy, 
Jo Jo 



/ / v^[A (* + 5(2/) - y) - F (t + 5 (y) - y)] dFi (y) - 
Jo JO 



/ dt U 1 (t + g{y)-y)dF 1 (y)= dF 1 (y) U^du 

Jo Jo Jo Jg(y) 



U 1 { U )F 1 (g- 1 (u))du, 




10 



V. Bagdonavicius; I. Masiulaityte and M. Nikulin 



We obtained 



where 



00 ft -J-. 

dt / F 1 (t + g(y)-y)d{VK[F 1 (y)-F 1 (y)]}-^ 
Jo 



POO l>t /"OO 

/ dt F 1 {t + g{y)-y)dU l {y) = F 2 (t)U 1 (t)dt- 
Jo Jo Jo 



poo poo 

I U 1 (y)[l-F 2 (y)]d(g(y)-y) = C/i (y) [rF 2 (y) - r + l]dy. 
Jo 



X°V 1 + V 2 + V 3 , 



f'OO POO 

Vi= / U(y)dy, V 2 = / %)£/i(y)dy, h(y) = c + r — 1 — F\(y/r) — rF 2 (y), 
Jo Jo 



000 

V3 = —rc 



POO 

/ U 2 (y)dy. 
Jo 



The variances of the random variables Vj are: 



Var(li) = Var(Ti), Var(F 3 ) = -^Var(T 2 ,) 

t2 



2 f°° f y 1 

Var(y 2 ) = - / [1 - F 1 (y)]h(y)dy / F 1 (z)h(z)dz = -Var(iJ(T u )), 
'1 jo jo '1 



where 



H{x) = / %)dy = x[c + r - 1 - F^x/r) - rF 2 (x)] + / ydF 1 (y/r)+ 
Jo Jo 



[ ydF 2 (y) =x[c + r-l- F^x/r) - rF 2 (x)) + 
Jo 



r 

lo 



rV(l {Tli < x/r} T Vi ) + rE(l {T2i < x} T 2i ). 

□ 

A consistent estimator of the variance a 2 is 

.. n ni ~ 2 ~ 2 "2 

- 2 = - - + 5 £[^) - ^ 2 + ^ £( T - - a 2 ) 2 , 

" i=l i=l n 2 i=1 
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where 



A = -VT 4 , c=— / 2/[l-F 2 (y)]dA(y) = - TT U [1 - F 2 (T U )], 

n ~l ^2 Jo ^2^1 ^ 



fl-(x) = x[c + f - 1 - A(x/f) - rF 2 (x)] + — J] l {Tl( <*/f}ri< + — Hn^x^i, 



ni <— ' ra 2 . 

i=i i=i 



«=i 

Theorem 3.2 : Suppose that rii/n — > Zj € (0, 1), n — > oo and £/ie densities fi(x), 
i = 1,2 are continuous and positive on (0, oo). Taen under Hq the statistic (14) 
converges in distribution to the normal law N(0,a 2 ), where 

a 2 = Var(Ti) + ^Var(H(T u )) + ivar(Q(T 2i )) 

'1 '2 

w/iere 

= Q(x)-^ 1 ( O - 1 (x)) +5 (x)[l-F 2 (x)]+E(l {9(Tii) < x} a(T 1 0)+E(l { ^<^ 

Q(x) = E{l {Tl( < x} [l - F 2 (T H )]//i(5(Ti,))}. 
Proof : Similarly as in Theorem 1 we obtain 

/ ^[F (1) (t) -F{t)]dt % \ U(t)dt, 
Jo Jo 

poo rt 

/ dt / y/K[F! (t + g(y) - y )-F 1 (t + g(y) - y)] dF l (y) - 

JO JO 

r vMy))-VM , , U1 „, ,,,, 

/"OO /'t 

/ / (t + g(y) -y)-F 1 (t + g(y) - y)\ dF l (y) - 

JO JO 



/>oo 

Jo 

/°° dt f F^t + g(y) - y)d{Mh{y) ~ F^y)]} Z 
Jo jo 
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/ dt F 1 (t + g(y)-y)dU 1 (y) = F 2 (t)U 1 (t)dt- 
Jo Jo Jo 



poo 

I7i(y)[i - F 2 (y)]d(5(y) - y) = / ^i(j/){^ 2 (y) - (</(y) - i)[i - ^(y)]}^. 

JO 







We obtained 



x°v 1 + v 2 + v 3 , 



where 

/'OO /"OC 

y 1= / U(y)dy, V 2 = / h{y)U 1 {y)dy, 
Jo Jo 



m = 7^)) i 1 - - ^(s-'to)) - *m») + (i/w - op - fid* 



The variances of the random variables V{ are: 

Var(Vi) = Var(ri), 



Va ,,n 2 [l-i^y)] 2 ^) / ^(^[l-i^]^^ 1 

Var(V3) = - / TTTT\ / 1 1 t \\ = r Var ^™)' 

h Jo /i(0U/)) Jo fi{9{z)) h 



2 r°° r y l 

Var(y 2 ) = - / [1 - F 1 (y)]h(y)dy / F 1 (z)h(z)dz = -Var(tf (T u )), 
n Jo Jo l i 

where 

Hlx) - £ ll uMr iFM ~ £ F ' (!r,fa)> * - £ 

f [i - WW*,) - [ [i - WW, = f (»)- 

F 1 {g-\x))x+ ydF 1 (g- 1 (y)) + [l-F 2 (x)]g(x)+ g(y)dF 2 (y) - x 
Jo Jo 

Jo [ \~My)) ] dFl[y) " + 9{x)[1 " F2(x)]+ 
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Sample size 


Significance level (%) 


50 


8.47 


100 


4.63 


170 


4.37 


200 


4.43 


400 


4.77 



n^MT^giTu)} + E[l {T2i < x} g(T 2i )} - x. 

□ 

A consistent estimator of the variance a 2 is 

-, n rai n 2 

- 2 = \ E< r « - ^ + S Eto) - ^ + J Eto) - <ft 2 > 

1=1 1 1=1 2 J=l 

where 

^ Hi 

= Q(x) - xAtrV)) + - F 2 (x)] + — J] l {§(Tl0 < ;r} 5(T li )+ 



^ n 2 ^ ni 

— E hT^xyg&i) - x, Q(x) = — i{T 1( <x}[l - ^(Tk)]//!^^)), 



r 1 (x)=F 2 " 1 (A(a:)), H = -Y j H{T u ), Q = -VQ(T 2iJ 
the density /i is estimated by the kernel estimator 



i=l 



4. Simulations 

We did a small simulation study of goodness-of-fit test for the hypothesis Hq. 
Failure times from exponential distribution were simulated: 

T y ~£(Ai), T 2j ~£(A 2 ), A 2 = rAi. 

The distribution function of a redundant system is 

F (t) = 1 - hl±^l e -^t + ^l e -(A 1+ A 2 )t_ 

Ai A 2 

The hypothesis //g is tested using 5 per cent significance level under several sample 
size n (see Table 1). Number of replications is 3000. 
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Tabic 2. Power of the test 



Sample size \ Constant 


0.1 


0.25 


0.5 


0.75 


100 


6 


12 


32 


85 


170 


11 


21 


79 


100 


400 


49 


88 


100 


100 



Alternative hypothesis Hq (at the switching time y the c.d.f. of the stand-by unit 
has a jump of size p(l — 

f { 2 y \x) = h{x + F^(F 2 (y) +p(l - F 2 (y)) - y) , < p < 1. 
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